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ABSTRACT 

We present the results from timing observations with the GMRT of the young pulsar 
J1833 — 1034, in the galactic supernova remnant G21.5— 0.9. We detect the presence of 4 
glitches in this pulsar over a period of 5.5 years, making it one of a set of pulsars that show 
fairly frequent glitches. The glitch amplitudes, characterized by the fractional change of the 
rotational frequency, range from 1 x 10~ 9 to 7 x 10~ 9 , with no evidence for any appreciable 
relaxation of the rotational frequency after the glitches. The fractional changes observed in 
the frequency derivative are of the order of 10 -5 . We show conclusively that, in spite of hav- 
ing significant timing noise, the sudden irregularities like glitches detected in this pulsar can 
not be modeled as smooth timing noise. Our timing solution also provides a stable estimate 
of the second derivative of the pulsar spin-down model, and a plausible value for the brak- 
ing index of 1.857, which, like the value for other such young pulsars, is much less than the 
canonical value of 3.0. PSR J1833 — 1034 appears to belong to a class of pulsars exhibiting 
fairly frequent occurrence of low amplitude glitches. This is further supported by an estimate 
of the glitch activity parameter, A g = 1.53 x 10~ 15 s~ 2 , which is found to be significantly 
lower than the trend of glitch activity versus characteristic age (or spin frequency derivative) 
that a majority of the glitching pulsars follow. We present evidence for a class of such young 
pulsars, including the Crab, where higher internal temperature of the neutron star could be 
responsible for the nature of the observed glitch activity. 
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1 INTRODUCTION 

Beside the basic, smooth spin-down of the neutron star due to the 
electromagnetic torque mechanism, pulsar timing studies also re- 
veal the presence of irregularities in the rotation of the star, mainly 
of two kinds : timing noise, which is characterized by continuous, 
random fluctuations in the rotation rate; and glitches, which are 
sudden increases in the rotation rate. These spin-up events are su- 
perposed on the long-term spin-down of the pulsar, and manifest 
themselves as sudden early arrival of the pulses. A recent work by 
fespinoza et al. ||201 ll) reports a total of 315 glitches observed in 
102 pulsars. The magnitude of the change in rotation frequency, v, 
during a glitch is typically in the range 1(T 10 < Av/i/ < 1(T 6 , 
and the fractional increment in the spin-down rate, Az>/i>, is in the 
range 1CT 5 to 10" 2 . 

The most plausible explanation for the sudden spin-up is the 
irregular flow of angular momentum from the faster rotating su- 
perfluid interior to the mor e slowly rotating s olid crust of the neu- 
tron star as it slows down jLvneetal.il 1995b . The current unified 
model for glitches is based on the superfluidity of the neutrons in 
a neutron star. The rotating superfluid in the neutron star carries 



angular momentum, by forming quantized vortices. The spacing 
between the vortices is negligible compared to the radius of the 
neutron star. On the macroscopic scale, the flow pattern looks like 
uniform rotation. These quantized vortices in the neutron superfluid 
in the inner crust can get pinned to the lattice of heavy neutron- 
rich nuclei. The pinning is possible because the effective width of 
the vortex core is less than or comparable to the lattice spacing of 
the nuclei. The pinning force is related to the energy gain when 
vortices are pinned to the lattice. The vortices stay pinned in this 
manner until a stronger force unpins them from the lattice sites. 
These pinned vortices in the crustal nuclei are rotating slower than 
the surrounding superfluid. Due to this differential velocity, mag- 
nus forces that act radially outward cause sudden unpinning and 
migration of vortices, which results in the transfer of angular mo- 
mentum from the superfluid to the crust. This gives rise to a sud- 
den speed-up of the solid crus t, which manifests as a gli tch in the 
timing behaviour of the pulsar. [Anderson and Itohl d 1975h were the 
first to make this connection between sudden unpinning of vortices 
and pulsar glitches. In the unpinned state, the superfluid moment of 
inertia is not coupled to the crust, hence the effective moment of in- 
ertia of the crust decreases, which in-turn increases the spin-down 
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rate. Between glitches, the vortex lines undergo a slow, thermally 
activated process, called vortex creep. The post-glitch relaxation is 
a process of recoupling of the vortices to another steady (pinned) 
state. Once the moment of inertia recovers due to this recoupling 
via repinning, the original extrapolated spin-down rate is restored. 
Thus, the observed sudden increase in the rotation rate, followed by 
exponential relaxation back to the extrapolated pre-glitch rotation 
rate, provides a useful probe of the neutron star interior. 

Th e pulsar J1833 — 103 4 was independe ntly discovered at the 
GMRT dGupta et al. I 120051) and at Parkes dCamilo et all l2006h 
and is associated with the galactic supernova remnant (SNR) 
G2 1.5— 0.9. This pulsar has quite a high spin-down luminosity that 
is amongst the top ten of all the known pulsars in our Galaxy. The 
flux density measured at radio wavelengths is very low — the es- 
timated mean flux density f rom the observations at 610 MHZ is 
0.65 mjy dGupta et al. Il2005h . With a period of 61.86 ms and a pe- 
riod deri vative of 2.02 x 10~ 13 s/s, it has a characteristic age, r c ~ 
4.8 kyr dCamilo et al7ll2006h , which makes it a fairly young pul- 
sar. Existing studies indicate that younger pulsars are more likely 
to show glitches. About half of all known pulsars with r c less than 
3x 10 4 have exhibite d glitches, but thi s fraction is much lower for 
the older population dYuan et al. l2010t) . PSR J1833- 1034 is thus a 
good candidate for the study of glitches. In this paper we report the 
detection of multiple glitches from this pulsar using timing obser- 
vations carried out at the GMRT at 610 MHz, and present a detailed 
study of its glitch activity. We also provide refined estimates of the 
timing parameters for this pulsar, including an estimate of the brak- 
ing index. In Section 2 we explain the observations and data anal- 
ysis techniques. Section 3 describes the detected glitches and their 
modeling in detail. In Section 4 we discuss the significance of our 
results. Summary and future scope are presented in section 5. 



2 OBSERVATIONS AND DATA ANALYSIS 

The GMRT is a multi-element aperture synthesis telescope con- 
sisting of 30 antennae, each of 45 m diame ter, spread over a re- 
gion of 25 km diameter dSwarupetaill997t) . Though designed to 
function primarily as an aperture synthesis telescope, the GMRT 
can also be used as an effective single dish in an array mode by 
adding the signals from indivi dual dishes, either coherently or in- 
coherently dGupta et al. I|2000T) . for studying compact objects like 
pulsars. The radio signals at the observing frequency band, from 
both polarizations of the 30 dishes, are eventually converted to 
baseband signals of 16 MHz or 32 MHz bandwidth, which are then 
sampled at Nyquist rate. These digitised signals are delay corrected 
and then Fourier Transformed in a FX correlator to get spectral in- 
formation. After fringe derotation, these dual polarization spectral 
voltage samples from all the antennae are added coherently in the 
GMRT Array Combiner (GAC), to produce the phased array out- 
puts for each polarization. These are then converted to intensity, 
integrated to the desired time constant and recorded on disk for off- 
line processing. The data are time-stamped using a minute pulse 
signal, derived from the observatory's GPS receiver, which is em- 
bedded in the data stream. 

The timing observations described here were carried out in the 
total intensity phased array mode at 610 MHz. In this mode of op- 
eration, the array needs to be phased up before observing the target 
pulsars. This is achieved by recording the correlator data for a point 
source calibrator, solving for the antenna based gains and phases 
from these, and applying the phases as corrections to the output of 
the Fourier Transform stage of the correlator. The array remains 



phased for up to a few hours and dephases due to slow changes in 
instrumental and ionospheric phases. When this happens, one needs 
to rephase the array to proceed with further observations. 

Timing observations for PSR J1833 — 1034 were started 
around mid-2005, shortly after its discovery at the GMRT. In the 
beginning, after the initial, closely spaced observations that are 
needed to obtain the timing solution for a newly discovered pul- 
sar, the observations were somewhat random and sparse in time. 
Since the occurrences of glitches are unpredictable and their relax- 
ation timescales can be quite short, regular monitoring is important 
for detection and study of the glitches. From mid-2007, after the 
possible detection of the first glitch from this pulsar, a regular tim- 
ing program was started, with observations roughly about 10 days 
apart, except for the GMRT maintenance intervals. Each observing 
epoch has a 90 min long scan on PSR J1833— 1034, and a shorter 
scan of 5 min on PSR B1855+09, which acts as a control pulsar 
for validating the data quality and reliability of the time-of-arrival 
(TOA) values from the newly established GMRT timing pipeline. 
The final data for each scan are total intensity values for each of 
256 spectral channels (across a 16 MHz bandwidth), recorded with 
a sampling interval of 0.256 ms. The main observing parameters 
for a typical epoch are summarized in Table[T] 

In the off-line processing, the recorded multi-channel total in- 
tensity data are first dedispersed to remove the effect of interstellar 
dispersion on the pulse shape. The dedispersed time series data are 
then synchronously folded using the topocentric pulsar period, ob- 
tained from the best existing model parameters (barycentric) for 
the concerned pulsar, after correcting for the observing time and 
location. The UTC corresponding to the middle of the observing 
session is used as the reference point for that particular epoch and 
it is derived from analysis of the GPS pulse signal embedded in 
the data. The predicted topocentric periods are calculated using 
"polyco" files produced by the pulsar timing program TEMPO q 
For the control pulsar, the barycentric model parameters were taken 
from the ATNF pulsar catalogQ; for J1833— 1034, these were ob- 
tained from the initial epochs of observations and refined at succes- 
sive epochs, as required. The topocentric TOAs at each epoch are 
obtained by cross-correlating the average profile at that epoch with 
the highest signal-to-noise (S/N) profile from all the epochs used 
as a template. These are then converted to solar system barycen- 
tric TOAs using the Jet Pro pulsion Laboratory DE200 solar system 
ephemeris ( Standishll982h inside TEMPO. To illustrate the quality 
of the profiles, Fig. Q] shows the highest S/N profile (which is used 
as the template), as well as a typical average profile (whose S/N 
is close the median value from profiles of all epochs). For the tim- 
ing analysis using TEMPO, we have used these topocentric TOAs 
along with the uncertainties related to the S/N of the profiles. The 
tiny uncertainty of the TOA of the reference epoch is artificially 
increased to make it close to the median value. 

At the solar system barycentre, the time evolution of the rota- 
tional phase of a solitary pulsar is well-approxi mated by a polyno- 
mial of the form dManchester and Taylor fl977h 

<M*) = 00 + v(t - to) + \v(t - to) 2 + -v{t - to) 3 (1) 

2 D 

where <f>o is the reference phase at time to; v, ii and v are the 
pulsar rotational frequency and its derivatives. TEMPO attempts 
to minimize the deviations between the observed and model rota- 
tional phases using \ 2 minimization. Timing irregularities are seen 

1 see http://www.atnf.csiro.au/researcli/pulsar/tempo 

2 see http://www.atnf.csiro.au/research/pulsar/psrcat 
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Table 1. Main parameters for the two pulsars observed, at 610 MHz with a bandwidth of 16 MHz, with a typical cadence of 10 days. 
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610MHz(mJy) 


time (min) 






B 1855+09 


5.36 


16.8+ 


5 


55970 


177 


J1833-1034 


61.86 


0.65* 


90 


87293 


32 



t extrapolated flux u sing the catalogued values at 400 and 1400 MHz. 
t lGupta et alH J2005I) 

N p is number of pulses accumulated in the integration time. 
(S/N) exp is the expected S/N for the 610 MHz profile at any epoch. 




Profile bins 

Figure 1. Highest signal-to-noise ratio profile (top panel) and typical av- 
erage profile (bottom panel) for PSR J 1833— 1034. The highest signal-to- 
noise ratio is 22, whereas the typical average profile has a signal-to-noise 
ratio of 10. 



as slow, large changes in the residuals (j> — 4>rn, where <j) is the 
measured phase and <j} m is the model phase. 

For a glitch, there is a sudden change in the phase residuals, 
modeled by an abrupt jump in the frequency and its derivative, fol- 
lowed by a relaxation proc ess. The frequency perturbation due to a 
glitch can be described as ( I Yuan et al71l2010l) 

Au(t) = Au p + AiJpt + Au d e~ t/Td (2) 

Au(t) = Avp + Au d e- t/Td (3) 

where Av and Ai> are the changes in the pulse frequency and its 
derivative, relative to the pre-glitch model; Av v and Av v are the 
persistent change in rotational frequency and its derivative; Av d 
is the amplitude of the exponentially decaying part of the jump in 
rotational frequency (and Av d is the corresponding value for the 
frequency derivative), with a relaxation time constant r d . The total 
frequency change at the time of the glitch is then given by 

Av 9 = Avp + Av d (4) 

The instantaneous change in v at the glitch is given by, 

Av g = Ai)p + Av d (5) 



3 RESULTS 

We first discuss the timing results for the control pulsar and then 
present the results from the timing analysis of the target pulsar 
J1833-1034, using GMRT data spanning 30th July 2005 to 11th 
Jan 201 1. The post-fit residuals obtained from the phase connected 
timing solution for the control pulsar B 1855+09 (shown in Fig. 
[2]l over the full data span of 52 epochs, yield a root-mean-square 
(rms) value of 15 ^is. This is more than the theoretical, expected 
value (ath) of 3 fj,s, which is based on the expected S/N of 177 
given in Table 1, and using <jtoa — -^j^ (where W is the pulse 
width and S/N is the expected signal-to-noise ratio of the average 
profiles). However, it matches well with the value expected for the 
achieved S/N (which has a maximum value of 35 in all the observ- 
ing epochs). Nevertheless, our final 610 MHz timing residuals for 
PSR B 1855+09 are wor se off by factor of 5 in rms from the 1420 
MHz results obtained bv lHobbs et all d2006h . In addition to signal 
to noise limitations, there can be effects of interstellar weather that 
can reduce the accuracy of the TOAs : the pulse arrival time at each 
epoch can have extra deviations due to propagation effects in the 
interstellar medium, which will be larger at the lower frequency. 
These results verify, to first order, the proper working of the pulsar 
timing set-up at the GMRT. 

In order to check for any low-level systematic effects in the 
timing residuals for this control pulsar, we investigated the changes 
in the rms value when adjacent residuals are averaged. For truly 
white-noise residuals, this rms should decrease as the square root 
of the number of residuals averaged. Fig. [3] shows the results for 
this on a log-log plot, where the data points are found to match 
quite well with a slope of —0.5 (green dashed line). The over-all 
post-fit residuals thus exhibit a white-noise behaviour and are likely 
free from any systematics. The averaging of post-fit residuals over 
166 days (10 TOAs) achieves a rms of 3 [is for B 1855+09, which 
implies a long-term timing stability of 1 part in 4x 10 12 . The above 
results for the control pulsar establish the basic fidelity of the timing 
pipeline for the GMRT. 

For PSR J1833-1034, since the S/N is typically significantly 
lower than for the control pulsar, to check the data quality for tim- 
ing purposes we show the distribution of TO A errors in Fig. [4] The 
bulk of the values are clustered in the range of 100 to 300 ^is, with 
a small tail of larger values. This skew in the distribution is due to 
degradation of S/N at some epochs, possibly due to fading caused 
by interstellar scintillation. However, the achieved TOA uncertain- 
ties are still good enough to detect changes in the residuals of the 
order of several milliseconds due to the occurrence of glitches. 

Starting with the initial timing observations for PSR 
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Figure 2. Timing residuals for the control MSP B 1855 +09, which show a 
rms ~ 15 fis. 
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Figure 5. Timing residuals for the target pulsar PSR J1833— 1034 from the 
first 1.5 years of data, prior to the first detected glitch at ~ 2007.2. 




Number of residuals averaged 

Figure 3. RMS of residuals versus number of consecutive residuals aver- 
aged, for the PSR B 1855+09. Green dashed lines indicate the expected 
slope of —0.5 for uncorrelated residuals. 
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Figure 4. Distribution of TOA errors for PSR J1833- 1034. 
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Figure 6. Residuals from the full 5.5 yrs of data without modeling of any 
glitch event, showing strong signature of timing noise, as well as evidence 
for glitches — these are typically seen as sudden negative change in the 
slope of the residuals, and the suspected locations are indicated by arrows. 
The detection and modeling of these glitches are explained in detail in the 
text, and illustrated in Figsr7ll8ll9l&[Tol 



J1833 — 1034, we are able to build up a phase-connected timing so- 
lution (shown in Fig. [5}, till the epoch of 2007.2. From this 1.5 
year data span, we obtain a fairly good timing model for this pul- 
sar, including a second frequency derivative (see the first row of 
Table |2), and the rms of the residuals is around 174 /xs. The ref- 
erence epoch (MJD) for these measurements is set to the epoch 
which is mid-point of our full data span (i.e. MJD of 54575), for 
better comparison with the later models. The pulsar position used 
in the ti ming model is the on e determined from the Chandra obser- 
vations dCamilo et al. l200fj) . The position derived from our timing 
solution of 1.5 years of phase-connected residuals is within the 3a 
error bars of this X-ray position. We derive a braking index (n — 
vvlii 1 ) of 2.168(8) for this pulsar from this initial data span (see 
last column of Table |2j. 

Fig. [6] shows the timing residuals for the full data set (94 
epochs spanning 5.5 years) for this pulsar, relative to a simple slow- 
down model including the pulsar spin frequency and its first two 
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Table 2. Rotational parameters for PSR J1833 — 1034 from different timing solutions. The first row is for the timing solution from the initial 1.5 years of data, 
before the first detected glitch. The second row is for the full 5.5 yrs of data, without inclusion of any glitch models. The third row is for the full 5.5 yrs of data, 
with glitch models fixed to the values derived from the piecemeal modeling of the glitches. The last row is for the solution from the final global fit, including 
all four glitches (2 free parameters each, with glitch epochs fixed to the values obtained from the piecemeal fittings), and four frequency derivatives. All errors 
are 1<t values. 
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22 
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2.168(8) 
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Figure 7. Timing residuals illustrating the first glitch event at ~ 2007.2. 
This initial model is obtained from fitting v, v and v to TOAs from 2005.5 
up to the suspected epoch of the glitch. The final model yields a glitch with 
Af s / v of 3.34 X 10 -9 , localised in time to the interval marked by the 
dotted lines. 
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Figure 8. Timing residuals illustrating the second glitch event at ~ 2007.9. 
This initial model is obtained from fitting v, i> and u and the glitch model for 
the first glitch, to TOAs from 2005.5 up to the suspected epoch of the glitch. 
The final model yields a glitch with Ai^ 9 / v of 1.00 X 10 -9 , localised in 
time to the interval marked by the dotted lines. 
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Figure 9. Timing residuals illustrating the third glitch event at ~ 2008.8. 
This initial model is obtained from fitting v & v over the first 130 days 
shown here, with v held constant at the value obtained from the fit to the first 
1.5 years of data. The final model yields a glitch with Au g / v of 1.6x 10 — 9 , 
localised in time to the interval marked by the dotted lines. 



Figure 10. Timing residuals illustrating the fourth glitch event at ~ 2009.9. 
This initial model is obtained from fitting v & v over the first 152 days 
shown here, with v held constant at the value obtained from the fit to the first 
1.5 years of data. The final model yields a glitch with Au g / v of 6.9 X 10 — 9 , 
localised in time to the interval marked by the dotted lines. 



derivatives. The best fit model parameters from this are listed in the 
second row of Table[2] The residuals, with a relatively large rms of 
15.4 ms, are clearly dominated by non-random, low frequency tim- 
ing noise effects. The amplitude of this timing noise is a strongly 
increasing function of the length of the data span. The effect of this 



timing noise was probably not detected for the initial data span of 
1.5 years, as the fitting of the spin-frequency and its two deriva- 
tives can mask most of the low-frequency trends. In these timing 
noise dominated residuals of Fig.[6j the presence of glitches can be 
distinguished by sudden changes in the slope of the curve. Clear 
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events are seen at 2007.2 and 2009.9, and less likely ones at 2007.9 
and 2008.8, all of which are marked by arrows in Fig. [6] 

The presence of a glitch is confirmed by taking relatively 
shorter stretches of data around the suspected glitch event, and do- 
ing a local timing fit to the TOAs, starting with a model for the data 
prior to the glitch. Sudden, systematic deviation of the residuals 
from a smooth behaviour is taken as the signature of the occur- 
rence of a glitch (as seen in Figs[71[8l[9l&|10ll. Detailed modeling 
is then carried out to estimate the glitch epoch, and the changes in 
frequency and frequency derivative at the glitch. The best possi- 
ble value of the glitch epoch is estimated by minimising the phase 
increment required to obtain a p hase-connected soluti on over the 
interval around the glitch epoch jjanssen et al. Il2006h . The mea- 
surement uncertainty of the glitch epoch is obtained from the cor- 
responding 3a limit of the glitch phase increment parameter. 

Starting with an initial timing model having v, v and i), we 
find that the first glitch (Fig. occurred at MJD = 51469 (± 7), 
with a fractional change in the rotational frequency (Avg/v) of 
3.34xl0" 9 . The modeling for this glitch includes 26 TOAs ob- 
served over 824 calender days, of which there are 236 days of 
data after the glitch event. Fig. [8] shows the second glitch event, 
which is best fit by a fractional increase in rotational frequency of 
l.OOx 10~ 9 at MJD = 54423 (± 9). As we have a fairly good tim- 
ing model for the first 1.5 years (including the first glitch event), 
without any timing noise effects, the pre-glitch interval for this sec- 
ond glitch includes all of these TOAs over 842 calender days, with a 
pre-fit model of v, v, v and the derived parameters of the first glitch. 
There is a third glitch (Fig.|9j detected at MJD = 54750 (± 15) with 
a fractional change in the rotational frequency of 1.6xl0 -9 . The 
last glitch event (Fig.QJjJ observed at MJD = 55142 (± 2) yields 
a fractional change in the rotational frequency of 6.9 x 10~ 9 . In or- 
der to minimse the effects of timing noise, the pre-glitch interval for 
the third glitch includes TOAs over 130 days and the fourth glitch 
includes TOAs over 152 days, with the pre-fit model having v, v 
& v. But since the pre-fit data span over smaller intervals, the fit 
uses v, ii as free parameters, with v kept constant to the value de- 
rived from the initial 1.5 yrs of data. Inclusion of TOAs over larger 
spans increases the influence of timing noise, where the residuals 
depart from the simple spin-down model with v, v and is, mak- 
ing it harder to detect the glitches accurately. There are also small 
changes in slow-down rate, of the order of 10 , observed at the 
glitch epochs. The new timing models, after inclusion of the glitch 
parameters, yield phase-connected timing residuals with rms val- 
ues of 177 fis, 216 ps, 227 /is and 1.5 ms respectively, for the four 
cases. 

It is sometimes possible that, for data that are relatively 
sparsely sampled and have significant amount of timing noise (both 
of which are somewhat true for the present case), there can be large 
deviations in residuals with respect to the basic spin-down model of 
v, v & v, which may mimic glitch-like behaviour. In order to dis- 
criminate between the effect of timing noise and genuine glitches, 
we investigated the pre-fit and post-fit residuals around the glitch 
epochs by fitting with higher frequency derivatives without inclu- 
sion of any glitch model. For example, for the case of TOAs span- 
ning over the first 936 days (including glitch- 1 and glitch-2), the 
rms for the post-fit residuals is 216 fis (shown in Fig, lilt. The 
model includes two glitches and a spin-down model with v, v & 
i>, which amount to 9 free parameters. The same span of TOAs can 
also be fitted with a model having v and 8 frequency derivatives, 
without inclusion of any glitch models, which also amounts to 9 
free parameters, as for the model with glitches. The post-fit residu- 
als shown in Fig.Q/2] have a much larger rms value of 951 fis and 
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Figure 11. Post-fit residuals after fitting for the first two glitches, for TOAs 
spanning the first 936 days. The rms is 216 fis. 
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Figure 12. Post-fit residuals (for TOAs spanning the first 936 days) after 
fitting with a model having v and 8 frequency derivatives, without inclusion 
of any glitch models, which amounts to 9 free parameters, as for the model 
with glitches. The rms is 95 1 /is, and there are large deviations at the epoch 
of the first glitch. 



also show large discontinuities, including at the glitch epoch. Sim- 
ilar effects were found for the data around the 3rd and 4th glitches. 
This illustrates the fact that some of the large TOA variations that 
we see for this pulsar, over and above the basic spin-down model, 
can not be satisfactorily explained with a model of timing noise 
characterised by higher order derivatives, but are better modeled 
with discrete glitch events. 

Finally, to obtain a global model for the entire data span, we 
have compared the following three approaches : (i) taking the mod- 
els for the 4 glitches obtained from the piecemeal fits to the individ- 
ual glitches as fixed and then fitting for v,v 8lv over the entire data 
span; (ii) a fit to the full data span using v and up to 12 frequency 
derivatives (the maximum allowed by TEMPO) without any glitch 
models included; and (iii) a global fit for 4 glitches, v and first four 
frequency derivatives, to achieve the same count of 13 free parame- 
ters as in case (ii) — each glitch contributes 2 free parameters (spin- 
frequency increment and change in spin-down rate), as the glitch 
epochs are fixed to the values obtained during the piecemeal fit- 
tings, by minimising the phase increment at the glitch epoch. Case 
(i) gives the residuals shown in the middle panel of Fig.[T3] with an 
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rms of 2.2 ms, and fairly smooth behaviour with large swings, typi- 
cal of timing noise. Results from this fit are given in the third row of 
Table |2] Case (ii) gives the residuals shown in the top panel of Fig. 
[T3l Though the rms is 1.4 ms, the variations of the residuals show 
sudden, large jumps (as at the epoch of the first glitch) and also 
sharp, cuspy variations (as at the epoch of the 4th glitch). For case 
(iii), global fits with all 4 glitches (using the results from the piece- 
meal fits as the starting pre-fit model) and increasing number of 
frequency derivatives were tried, and the following was found : for 
4 glitches plus v,v the residuals are 1.2 ms and the behaviour 
is qualitatively similar to case (i); for the case of two more deriva- 
tives added to achieve the same number of 13 degrees of freedom 
as it was for the case (ii), the residuals (shown in the bottom panel 
of Fig. 113b reduce significantly to 0.5 ms and also the slow, large 
fluctuations typical of timing noise are noticeably suppressed. The 
glitch parameters are not very different from those obtained from 
the piecemeal fits. The results from this model are summarised in 
the 4th row of Tableland the final glitch parameters are given in 
Tabled 

From the above, we argue that the best timing model is that 
given by a global fit of 4 glitches and 5 spin frequency terms, which 
gives the best global fit to the data and reduces the residuals to a 
minimum. The attempt to fit the TOAs with a pure timing noise 
model having large number of derivatives does not give accept- 
able results : both for localised fits to data sets that span individual 
glitches, as well as for the global data set. For such cases, the rms 
of the residuals is larger and/or the residuals show uncharacteristi- 
cally large, abrupt changes. We take the results from this fit as the 
final timing model for this data set. These results are summarised 
in the last row of Tableland in Table[3] 

Now in order to measure the amount of timing noise 
present in this pulsar, w e have used the definition given by 
lArzoumanian et al71 dl994l) . 



A 8 = 



(6) 



where the spin-frequency, v and its second derivatives, v, are mea- 
sured over at — 10 8 s interval. We have used first 3.16 years of data 
for PSR J1833- 1034 to estimate the value of A 8 as 0.5, which fol- 
lows the correlation between timing noise and spin-down rate, i.e. 
the younger pulsars with larger spi n-down rate exhibit more ti ming 
noise t han older pulsars, se en by lArzoumanian et al71 i 19941) and 
later bv lHobbs et aTI feoid) . 

The value of the braking index determined from the final 
global fit is 1.8569(6). This braking index is much less than 3, 
which is in general agreement with the values obtained for other 
young pulsars having reliable estimates for this quantity. For ex - 
ample, for Crab pulsar, n = 2.509(1) dLvne et al. 1119881 ^~ 
for PSR J1846-0258, n = 2.65( 1) (Livingstone et aT 



1993) 



2006), 



for PSR B0540-69, n = 2.140( 9) ^Livingstone et al.ll2005l) . for 
PSR B 1509-58, n = 2.837 (1) dKaspi et al. Ill994h and for PSR 
Jl 119-6127, n = 2.684(2) dWeltevrede et al. 11201 ll) . A value of 



n < 3 indicates that simple magnetic dipole model does not com- 
pletely explain spin-down evolution of pulsars. Particle outflow in 
the pulsar wind can also carry away some of its rotational kinetic 
energy. 



4 DISCUSSION 

Our timing study of the young pulsar J1833 — 1034 associated with 
the galactic SNR G21.5— 0.9 shows clear evidence of frequent 
glitches in the pulsar's rotational history. We find as many as 4 
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Figure 13. Timing residuals for PSR J1833— 1034, for different types of 
global fits, all having the same number of 13 free parameters. The top panel 
shows residuals with 12 frequency derivatives fitted, without any glitch 
models. The suspected glitch epochs are marked by arrows. At the first, third 
and fourth glitch epochs there is clear evidence for discontinuities in the 
residuals. The middle panel shows residuals with a model of four glitches, 
derived from piecemeal fitting and two derivatives of pulsar spin-frequency. 
The presence of low-frequency timing noise is clearly seen in these resid- 
uals, which otherwise show a smooth behaviour. The bottom panel shows 
residuals obtained from a global fit of four glitch models and four frequency 
derivatives. The rms is significantly lower and the residuals are much more 
whitened in nature, indicating that this model gives the best global timing 
fit to this data set. 
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Table 3. The parameters of all the four glitches detected in PSR 1833—1034, as determined from a global fit to the timing data. The errors on the least 
significant digit are at la level. 



Glitch epoch 
(MJD) 


Date 


Fit span 
(MJD) 


AVg 

(10- 9 ) 


Ai> a 

(io- 5 ) 


54169 (±7) 


5th Mar2007 


53581-54405 


3.11(5) 


1.4(2) 


54423 (±9) 


12th Nov2007 


53581-54517 


1.09(6) 


4.0(3) 


54750 (±15) 


1 1th Nov2008 


54620-54885 


3.55(6) 


-7.7(2) 


55142 (±2) 


6th Nov2009 


54990-55572 


7.50(8) 


-9.9(2) 



] 0.001 

Change 



0.01 0.1 1 10 100 1000 

in spin-down rate (1E-15 Hz s A — 1) 



Figure 14. Plot of changes in the pulsar spin-down rate, \Ai>\, with spin- 
frequency due to glitches. The data points (red circles) are from the glitch 
table of the ATNF pulsar catalog. The data points for PSR 1833-1034 are 
denoted by squares. 



glitches over the observing span of 5.5 years. Compared to the typi- 
cal range of glitch amplitudes mentioned in section 1, the fractional 
changes in the rotational frequency seen for this pulsar are rela- 
tively small, ranging from 1 x 1CP 9 to 7 x 1CP 9 . This behaviour 
is similar to the Crab pulsar, which shows Av/v ~ 10~ 8 ; whereas 
the Vela pulsar exhibits larger glitches, generally with Av/v > 
1CP 6 . As the amplitude of a glitch is related to the amount of stress 
built up in the pinned vortices, one might expect some correlation 
between the amplitude of glitches and the inter-glitch interval. Pul- 
sars that have small amplitude glitches do tend to show smaller 
interval between glitches (as observed for PSR J053 7— 6910 by 
iMiddleditch et al. 1 J2006h and for PSR B 1642 -03 by IShabanoval 
1 20091) ). and this is borne out in the case of PSR J1833 — 1034 as 
well. Clearly, this pulsar falls under the category of pulsars that 
exhibit relatively frequent, but low amplitude glitches. 

PSR J1833 — 1034 also shows small but permanent changes in 
the slow-down rate at the glitches, and the typical fractional change 



of v is a few parts in 10 (Table[3]l. These small increases in 
are thought to be due to the decrease in the effective moment of in- 
ertia of the crust, which includes all components of the star dynam- 
ically coupled to the crust. Decoupling of the superfluid moment 
of inertia during the unpinning state reduces the entire moment of 
inertia of the star. However, for the third and fourth glitches, we ob- 
serve a decrease in v. This sign change of v may imply a small in- 
crease in moment of inertia or a small decrease in spin-down torque 
at the time of the glitch. 

We did not detect any exponential recovery or decay with time 
after the glitches in our data, for either the change in rotational fre- 
quency or its derivative. This may imply that there are only perma- 



nent changes in the rotational parameters when this pulsar glitches. 
However, there is also a possibility that this may be due to the fact 
that our sampling interval for the timing properties of this pulsar 
— about a week to 10 days — is somewhat coarser than what may 
be required to adequately sample the expected decay time-scale for 
such low amplitude glitches. For example, in case of the Crab pul- 
sar, for glitches with an amplitude of the o rder of ~ 10~ s , the expo- 
nential decay time-scale is of ~ 10 days dWong etal.ll200ll) . Such 
time scales would be hard to detect in our timing data, and would 
need a much more intensive campaign of observations. 

For the general pulsar population it is found that glitches with 
small Av also have small changes in \i>\. This is shown in Fig. 1 141 
using the database of the glitch table in the ATNF pulsar catalog q 
where a clear correlated trend can be seen. Our results of the glitch 
parameters for PSR J 1833— 1034 shows that this pulsar follows this 
trend quite well. 

The level of strength and frequency of occurrence of glitches 
in a pulsar can be quantified by the glitch activity parameter, A g , 
defined as the mean change in frequency per unit time owing to 
glitches dLvne II 19991) : 



(7) 



where ^ v g ' s the total increase of the frequency owing to all the 
glitches over an interval of T. Glitches are considered as events of 
angular momentum transfer from the superfluid interior to the crust 
of the neutron star. The same rate of angular momentum transfer 
can be achieved with frequent small glitches or occasional larger 
ones. The glitch activity parameter combines the amplitude and 
frequency of angular momentum loss due to glitches over the in- 
terval of T. A g is relatively insensitive to the additional discov- 
ery of smaller glitches as the quality of a given data set improves, 
and hence it can be used as a long-term indicator of glitch effects 
dWong etal.ll200lh . We find A g = 1.53 x 10 -15 s~ 2 for PSR 
J1833-1034. 

Fig. [T5] shows the range of known values of A g , as well 
as its dependence on \i>\, for a collection of 32 pulsars. The 
data are mostly from literature (circles), except for a few points 
(triangles for B0611+22, B 1853+01 and B0540-69) which are 
from unpublished results from observations at the Torun Radio 
Telescope by one of us (Wojc iech Lewandowski ). The litera- 
ture references are as follow : iLvne et al. I d2000|) (B0833-45, 
B1325-43, B1535-56, B1641-45, B1727-33, B1736-29, 
B1758-23, B1800-21, B1823-13, B1830-08, B 1859+07 
B2224 +65, B0355+54, B0525+21 and B 1737-30), IWang et al. I 
J2000h (B0833-45, B1046-58, Jl 105-6107, Jl 123-6259, 
B1338-62, B1610-50, B1706-44, B1727-47, B1758-23, 



3 see http://www.atnf.csiro.au/research/pulsar/psrcat/glitchTbl.html 
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Figure 15. Plot of glitch activity parameter, A g , versus the spin-down rate, 
i>, of pulsars. Circles are data points taken from literature and triangles 
are data points taken from Torun observations. The glitch activity of PSR 
1833—1034 is denoted by the square. 



Characteristic age 
Figure 16. Plot of glitch activity parameter, A g , 



(kyr) 

versus characteristic age 



of pulsars. Circles are data points taken from literature and triangles are 
data points taken from Torun observations. The glitch activity of PSR 
1833—1034 is denoted by the square. 



B1757-24 and B1 800-21,), IWong et al. N200ll) (B0531+21), 
Hobbs etal.N2002t) (J18 06-2125). lUramal d2002l) (B1737-30) 
Weltevrede et al. I J201 lh (Jl 119-6127), [Shabanova andUramal 



fcOOOh (B1822-09). iMiddleditch et al~l d2006l) ) (J0537-6910), 
iLivingstone et al. I d2006l) (J1846-0258). Though there is some 
scatter present, for a majority of pulsars there is an overall trend of 
increasing A g with increasing This trend is mirrored in a plot of 
glitch activity versus characteristic age, as shown in Fig.[T6]: glitch 
activity is higher for younger pulsars with characteristic age ~ 10 
kyr, and as the characteristic age increases, the activity falls off. 
These effects could be due to the fact that the flow of the angular 
momentum from the interior decreases with age (or increases with 
|i>|). There are a few pulsars with relatively higher values of \v\ (or 
relatively smaller values of characteristic age) that have somewhat 
lower values of A g , and hence lie off the main curve. Detailed in- 
vestigation shows that these are a group of very young pulsars (i.e. 
low characteristic age), such as the Crab, PSR Jl 119— 6127, PSR 
B 1853+01, PSR J1846-0258 and PSR B0540-69. We find that 
our young pulsar J1833— 1034 fits in very well with this group. 

Glitches are thought to be caused by the release of stress built 
up during the regular spin-down of the pulsar. This stress on the 
pinned vortices in the superfluid interior gets released to the solid 
crust by a collective unpinning of many vortices. This unpinning 
process results in a sudden spin-up of the crust due to this discontin- 
uous transfer of angular momentum from the interior, which in-turn 
is manifested in a change in observed pulsar frequency. Frequent, 
low amplitude glitches implies that the release of the built up stress 
happens in a more uniform and continuous manner than for pulsars 
which show few, large amplitude glitches. In other words, for the 
younger pulsars with larger slow-down rates, the flow of the angu- 
lar momentu m from the interior seems t o be a smoother process. 
According to lMcKenna and Lvne I ( fl990l) . the higher internal tem- 
perature associated with the younger neutron stars might prevent 
the build up of larger stresses. In such cases, stresses on the pinned 
vortices get relieved by thermal drift of the vortices from one pin- 
ning site to another in a gradual fashion, resulting in frequent low 
amplitude glitches. Hence such pulsars may constitute a distinct 
class of glitching pulsars: younger pulsars with lower glitch activ- 
ity and higher internal temperatures. These relatively young pulsars 
may evolve towards the normal trend (i.e. towards right on Fig. [16} 
as they age. 



5 SUMMARY AND FUTURE SCOPE 

In this paper we have presented results for four glitches detected 
in PSR J1833 — 1034, from 5.5 years of timing observations at the 
GMRT. These glitches show fractional change of the rotational fre- 
quency ranging from 1 x 10 -9 to 7 x 10 -9 , with no evidence for any 
appreciable relaxation of the rotational frequency after the glitches. 
The fractional changes observed in the frequency derivative for this 
pulsar are of the order of 10 . This pulsar appears to belong to a 
class of pulsars exhibiting fairly frequent occurrences of low am- 
plitude glitches. We calculate the glitch activity parameter for PSR 
J1833-1034 to be 1.53x 10~ 15 s~ 2 , which puts it in a special class 
of young pulsars like the Crab, and offset from the normal trend of 
glitch activity versus characteristic age (or spin frequency deriva- 
tive) that a majority of the glitching pulsars follow. This could be 
related to the thermal history of young neutron stars. 

The final timing solution obtained after modeling of the 
glitches provides reliable estimates of the second derivative of the 
spin-down model for PSR J1833 — 1034. The resulting braking in- 
dex of 1.8569(6) is much less than the canonical value of 3, as also 
found for other young pulsars, supports the claim that pure dipole 
braking does not provide the full picture for pulsar spin-down. 

With aid of the high time resolution and coheren t dedisper- 
sion c apabilities of the new GMRT Software Backend l lRov et alTl 
l2010l) . we aim to search of giant pulse (GP) emission from this 
young pulsar. Though it is thought that most of the GP emitters 
are neutron stars with stro ng magnetic field at the lig ht cylin- 
der (B LC = 10 4 to 10 5 G) teomani and Johnston I l200lf). the de- 
tectio n of GPs in p ulsars like J 1752+2359 jErshov and Kuzminl 
120061), Bll 12+50 tershov and KuzmTnl J2003l) and B0031-07 
dKuzmin et al. 1120041 ; iKuzmin and Ershov 112004b reveals that GPs 
are also produced in pulsars with relatively low magnetic fields at 
the light cylinder. So even though the Blc of J1833-1034 is factor 
of 7 lower than the Crab, it can be worth searching for GPs using 
the coherent dedispersed output taken with the GSB. 
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